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We formulate the problem of a two-level system in a linearly polarized laser field 
in terms of a nonlinear Riccati-type differential equation and solve the equation an- 
alytically in time intervals much shorter than half the optical period. The analytical 
'sj" I solutions for subsequent intervals are then stuck together in an iterative procedure 

^ ' to cover the scale time of the laser pulse. This approach is applicable to pulses of 

^SJ ' arbitrary (nonrelativistic) strengths, shapes and durations, thus covering the whole 

region of light-matter couplings from weak through moderate to strong ones. The 
method allows quick insight into different problems from the field of light-matter 
00 i interaction. Very good quality of the method is shown by recovering with it a num- 

ber of subtle effects met in earlier numerically calculated photon-emission spectra 
from model molecular ions, double quantum wells, atoms and semiconductors. The 
method presented is an efficient mathematical tool to describe novel effects in the 
region of, e.g., extreme nonlinear optics, i.e., when two-level systems are exposed to 
pulses of only a few cycles in duration and strength ensuring the Rabi frequency to 
^ ' approach and even exceed the laser light frequence. 

^ ' 

O ; PACS numbers: 42.50.Md, 42.50.Hz, 42.65.Ky, 42.65.Re 

I. INTRODUCTION 

^ , 

^ , In the theory of light-matter interactions there is probably no more fundamental model 

^1 than the two-level one j^. Over the last decade, for example, the model has succedded 
^ ^ in explaining the main features of propagation of strong a few-cycle pulses through atomic 

^ . and semiconductor media 0, 0, 0, IE ific i S HO, 11, 113 ; e.g., carrier wave Rabi flop- 



^ ! ping, third-harmonic generation in disguise of second harmonic and carrier-envelope phase 
effects, to name a few. It has also made a basis for the descri ptio n of high-order harmonic 
generation from a single ato m jT3 | . a symmetric molecular ion |l4l . isl 3| and a double- well 
quantum structure 17, [3, Il9l| with emphasis on the strongly non-perturbative picture of 
the phenomenon and the occurrence of peaks in the spectrum of coherently scattered light 
at the positions of even harmonics. When applied to double wells, the model turned out to 
be successful also in explaining the effect of laser control of tunneling [2o| . 

Despite its dissemination in atomic, molecular and solid state physics the two-level model 
of light-matter interaction still suffers from the lack of exact analytical solution covering the 
whole range of laser intensities as well as pulse shapes and durations. The analytical solutions 
known hitherto cover only some different limiting cases. For instance, the most celebrated 
rotating-wave-approximation (RWA) solution is restricted to laser intensities ensuring 
the resonance Rabi frequency, Qr, to be much smaller than the laser frequency, u. Beyond 
RWA, the known analytical solutions include the non-RWA corrections along a perturbative 
procedure (e.g. j2l|, H^), either are valid in the so-called multiphoton excitation region 
[Til l {u « ijj2\ along with « U21, where UJ21 stands for the frequency separation 
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II THE ANALITYCAL SOLUTION 



between the two levels) or in the quite opposite strong coupling region ij, IE 3, 12 , 3, 1^ 
(a; >> LJ21 and Qr » u!2i)- Probably, the only analytical solution covering the whole 
intensity region is the recent one of Tritschler, Miicke and Wegener jUj for a box-shaped 
pulse, but obtained within the so-called square-wave approximation (SWA) consisting in 
replacing the actual time behavior of the field within half the optical cycle by a square of 
a constant appropriately chosen magnitude. Being approximate, this solution was able to 
reproduce only qualitatively some features of the exact numerical calculations, especially 
for the case of the resonant excitation (a; = U21), but was less convincing for distinctly 
off-resonant excitation. 

The aim of our paper is to present a quick iterative procedure for the problem of the 
two-level system in linearly polarized laser field, based on an analytical solution of the 
Schrodinger equation in very short time intervals. The analytical solution turned out to be 
possible thanks to defining the problem of level populations in terms of a single nonlinear 
Riccati-type differential equation in conjunction with dividing each halfcycle of the pulse 
into a number of narrow slices of equal width and considering constant the electric field 
within each slice with a value determined by the pulse function at the middle of the slice. 
The analytically obtained solutions for all slices in the pulse are then stick together by a 
simple recurrence formula derived relating the boundary conditions in the adjacent slices. 
This approach offers a simple analytical formula for the ratio of level population amplitudes 
within each slice, resulting in equally simple analytical formulae for population inversion, 
induced dipole moment and spectrum of the radiation emitted by this dipole. The photon- 
emission spectra obtained al ong the above line reproduce the numerically calculated ones 
available in literature [1, ^ lisi Eil . 17, [3, [3- Moreover , o ur iterative method 

indicates weak points of the square- wave solution of Tritschler at el. [23[ and is proved to be 
particularly useful in the area of extreme nonlinear optics |23||], i.e., when a few-cycle, strong 
pulses stimulate significant population dynamics in a two-level system on a time scale of 
half the optical cycle. 



II. THE ANALYTICAL SOLUTION FOR SHORT TIME INTERVALS AND 

ITERATIVE METHOD 

When presenting our analytical solution for short time intervals we start with the standard 
expansion il^lt) = bi{t)\l) + 62(^)12) for the wave function of the two- level system in a laser 
field, where |1) and |2) stands for the time-independent opposite-parity eigenstates of the 
bare system with eigenfrequencies Ui and U2, respectively. The time-dependent population 
amplitudes of the levels, bi{t) and 62(^)5 are then governed [l^ by the equation 

i-^bk = ujkh - ^it)bi{t), (1) 

where both k and / run the values 1, 2 with the constrain / 7^ k, and Q(t) = Q[ih{t) is the 
instantaneous Rabi frequency with Qji = fieo/h being the usual Rabi frequency as deter- 
mined by the dipole transition matrix element /i = (l|e2;|2) and the electric field amplitude 
eo, while h(t) = f(t)sin{u!t + 0) describes the incident-field evolution with f(t) having the 
sense of pulse shape (for pulses of at least few cycles in duration), a; the carrier frequency 
and the carrier-envelope offset phase. The latter is known [1, [l^, 0, H^l to be a 

relevant quantity determining the response of the system in the regime of few-cycle pulses. 
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Traditionally, it has been solved in different coupling regirnes either a set of two linear 
differential equations for bk with no RWA applied (e.g Iisl EgL U^) or more often (e.g. 
jisi E3 . Esl . tl9„ .23.] ) the resulting set of three linear differential equations for the Bloch 
vector components: u = 2Re(b^b2), v = 2Im(b^b2), w = — l^iP- Instead, we prefer 
to work with only one but nonlinear differential equation for the ratio r{t) = b2{t)/bi(t) of 
the population amplitudes. Through the population conservation law, + = 1, the 
above r determines directly both the population inversion w = (|rp - l)/(|r|2 + 1) and the 
induced dipole moment d{t) = {ip{t)\ez\ip{t)) = fiu = 2/iRe(r)/(|rp + 1) and, consequently, 
the spectrum of coherently scattered light as well. After differentiating r over time and then 
using Eq. (P) one obtains [2l[ |2^ r to fulfil the following differential equation: 

i^ = ir^-lMt)+^2ir, (2) 

where U21 = uj2 — uji is the frequency separation between the bare levels. This equation 
falls into the family of nonlinear Riccati-type equations and a way for its iterative solution 
results from the transformation 

r(t) = ^^{n^ff{t)Rit)-^2i). (3) 



n^ff{t) = JAn\t) + ui, (4) 



converting Eq. Q into 



/T52 ^\^^^^ . ■('^21 r-, 1^ ^21 



To avoid the cumbersome second term on the right-hand side, including the derivative 
dQ/dt, we divide the time scale of the pulse into a number of sufficiently narrow intervals 
with fj < tj < tj being the running time within the jth interval. In each interval of its 
width much shorter than half an optical cycle we approximate the Rabi frequencies ^{tj) 
and VL^^^ (tj) as constants of the values which they actually take in the middle (t™) of the 
interval. Under such an approximation, Eq. (0) when adapted to the jth interval looks 
hke i{dRj/dtj) = [R] - l)nf^/2, where Vlf^ = Q^ff{tj = tf). The resulting equation has 
straightforward analytical solution 



1 cot (n'/^tj -t))/2) 

R,{tj) = (6) 

Rf-i cot [nf^{tj~ti)/2j 



where i?*" = Rj{tj = t*) is the initial value of Rj, i.e., that at the beginning of the jth 

t-j, we have Rjitj = t-j) = -R*" 1, 



interval. For the extreme time in the interval, t, = t^ we have = t^A = R^,, resulting 



in the recurrence formula for the initial conditions 



1 -iRf cot (n'/^{t^^ -ti)/2) 

K+i = — ^ — - — (7) 

Rf-tcot (^fif^(tj-t;.)/2j 
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As a consequence of equations © and ((Tj) we obtain from Eq. Q the solution for rj(tj] 

2Qj - (u2i + inf^ cot (nf^itj - t))/2 



UJ21 



iQlj" cot (nf^{tj - t])/2) + 2fi 



(8) 



and also the recurrence formula for the initial conditions, r*", at the beginnings of subsequent 
time intervals 



2^- - (uj2i + il^f ^ cot - t;)/2)) 

U21 - inf cot fof ^(tj - t])/2) + 2^r;" 



(9) 



where Vf^^^ 



An] 



UJ21 with Qj = QrHj and /ij = f{t™') sin (cut™' +0). The solutions in the 

form of Eqs (jH)) and © allow us to obtain population inversion, induced dipole moment and 
electric field of coherently scattered light within the subsequent time intervals, < tj < tj, 
and to stick the solutions for the intervals to cover the whole time scale of the incident pulse. 

Before writing down the final formulae it is convenient to introduce the dimensionless 
strength parameter x = Qpi/u}, the dimensionless level separation parameter y = UJ21/LJ 
and the dimensionless time parameter r = ut, where < r < 2ttN for a A^-cycle pulse. 
Then, we divide each halfcycle in the r domain into K intervals of width n/K each, letting 
j to fall into the range 1 < j < 2NK. Within the jth interval, the running time covers 
the range {j — = rj < tj < rj = jj^ and the middle of the interval occurs at r™ = 
^T< ~ 27? • ^Iso, we make the replacement r*" = Ij and introduce the normalized effective 

^ jijj = ^JAx] + y^, where Xj = xhj with 



Rabi frequency within the jth interval 



eff 



eff 



f) sin(rf + 



) . In this language the recurrence formula of Eq. (jH)) reads 



y + ix. 



cot (^7rxfV2^)) ^: 



f ^ cot ( TTxf V2i^ ) + 2Xj/j 



ix. 



(10) 



For a given field-system parameters x,y,(j), /{r), Eq. allows us to generate the initial 
conditions for all subsequent 2NK time intervals from the only known initial condition Ii 
for the first interval (Ji = throughout this paper). Having generated the initial conditions 
we calculate the evolution of population inversion within the j'th time interval from 



y 



+ 4xjRe(Ij)) 



+4x, (x,(l - |/,f ) - yRe(Ij)) cos(xf ^(r, - rj)) - Ax.xflv^il,) sin(xf (r^ - , (11) 



while the evolution of the induced dipole moment from 

d-ir ) = X 

/I I \T \2\( eff\2 



(2/(l-|/,f)+4x,Re(Ij)) 



-y (x,(l - |J,f ) - yRe{l,)) cos(xf ^(r, - rj)) + yxflm{l,) sin(xf (r^ - ij))] , (12) 
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where < — rj < n/K within each intervaL Taking the second derivative of Eq. ()12|) 
with respect to Tj results in the electric field of the coherently scattered light in the forward 
direction: 



l + |/,f 



X 



yReil,)) cos{xf{T, - r])) - xflm{l,) sin(xf (r^ - rj)) 



(13) 



To study spectra we need to take Fourier transforms {tj z, where z is the spectrometer 
frequency in units of the incident light frequency u) of equations (fT^ and (fT^ with the 
results 



dj{z) 



fie 



-iz{j — 1)tt/ K 



|J,f)+4x,Re(/,))2/o 



-y {x,{l - |/,f ) - yRe(/,)) (/ri + //i) ) + yxf Im(/,)(/ri - / 



(14) 



and 



fiye 



-iz{j — l)iT/ K 



l + |/,f 



where 



(x,(l - |/,n - i/Re(/,)) (/ri + /;i) - xf Im(/,)(/ri - // 

^-i{z+qx''.^^)TT/K _ 2 



I e// 

2; + gXj- 



(15) 



(16) 



with q = 0, ±1. Finally, to cover the whole time scale of the pulse one needs to sum up 
equations (fTT |) - (fT3j) over j, taking into account equations (fTUj) and (fT^ . 



III. QUALITY OF THE ITERATIVE METHOD 

We have extensively examined the accuracy of the iterative method (equations ()lUp ~ 
(fTBj) ) in wide ranges of pulse shapes /(t), pulse strengths x, carrier frequencies y and carrier- 
envelope phases 0. In any case the method was found to be able to fit the results of direct 
numerical integrations of the Riccati-type Eq. Q, provided that K, i.e., the number of 
intervals into which we divide each optical halfcycle was chosen appropriately. Generally, 
the higher K the better was the quality of the method, as expected. However, K of the order 
of only a few units or at most ten appeared to be sufficient to ensure good-quality of the 
method for not too strong pulses (x < 1). For extremely strong pulses {x >> 1), generating 
fast population dynamics on a time scale of an optical cycle, an increase in K was needed for 
the method to reproduce all details of the numerical solution. However, even in the latter 
case only a little of computer time was consumed to accomplish successfully the iterative 
procedure with the use of the short-time-interval analytical solutions, i.e., equations (jlUj) - 

(HSl). 

To exemplify the effect of better quality of the iterative method when increasing K, let 
us focus on the one-photon resonance {y = 1) by a pulse of moderate strength (x = 1). 
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III QUALITY OF THE ITERATIVE METHOD 



We intentionally take this case because it is covered neither by the strong-coupling {y « 1 
and X » y) analytical solution of Ivanov et al. jlil . 17 1 nor by the multiphoton-excitation 



Dv tne n 

{y » 1 and x « y) analytical solution of Zuo et al. [3 [13 • Moreover, to assess 
the square-wave-approximation (SWA) solution of Tritschler et al. [i^ we choose a box- 
shaped (/(r) = 1) sine-like (0 = 0) pulse. The SWA, originally applied to the system 
of optical Bloch equations, consisted in replacing the sequence of halfcycles of the electric 
field by the sequence of identical squares, each of a height ensuring the areas under the 
halfcycle and square to be equal. In terms of our short-time-interval analytical solution, 
SWA corresponds to the choice of = 1 {hj = (—1)-'^^) and to rescaling x —>■ ^x resulting 
in Xj = {—ly^'^^x. In this limit our equations for Wjijj) and djijj) convert into those of 
Tritschler et al. obtained by a different analytical approach exploiting the Bloch equations. 
For the pulse of = 2 cycles in duration, now available in the laboratory practice (e.g. (23^), 
we show in Fig. (^) the effect of K on the population inversion calculated iteratively with 
the use of Eq. (fTTj) (solid lines), and compare this result with that obtained by integrating 
numerically the Riccati-type Eq. Q for r (dotted line). As seen, the choice oi K = 
ensures nearly perfect coincidence between the two approaches. On the other hand. Fig. 
(HJd) provides a comparison between our iterative results at iiT = 10 and the SWA results 
(dashed line) leading to a conclusion that the square-wave approximation can be used only 
to general predictions of qualitative nature. 

To prove a good quality of our iterative method we now recover with it some numerically 
calculated spectra of light coherently scattered by two-level systems, available in the litera- 
ture. One such a two-level system that has received a lot of attention in the past is the lowest 
pair of different-symmetry electronic levels of the H2 molecular ion (lo"g and Icr^), a pair be- 
ing well isolated from other levels particularly at large internuclear distances. In particular, 
Zuo et al. show in Fig. 6b of their paper 3| the two-level numerically calculated photon- 
emission spectrum of H2 in the near-resonance region translating into our y = 1.1 and 
X = 1.86. The spectrum was obtained by assuming the /(r) cos(r) electric field with /(r) 
gaussian increasing by 10 optical cycles (/(t) = exp[— ((r— 207r)/107r)^] for < r < 207r) and 
then keeping a constant value up to 30 cycles (/(r) = 1 for 207r < r < 607r). For the above 
conditions, we apply our Eq. (fT^ (with put to 1) along with Eq. (fTUj) to present in Fig. 
^) the iteratively calculated spectrum \d{z)\'^ = \ Y.jdj{z)\'^ with 1 < 3 < 2NK = 60K. 
To achieve high resolution of our spectrum we chose K = 100 and we will maintain this 
choice through all other figures to be presented. Our spectrum of Fig. (j2K) consists of 
Mollow triplets occurring at each odd-order harmonic (1, 3, 5 and 7) with the same sideband 
separation within the triplets. This iteratively obtained structure is in full agreement with 
the numerical spectrum of Zuo et al. (Fig. 6b in |1^]). In a different paper Zuo et al. 3 
give, for the /(r) cost = cosr field, their numerical spectrum for the same system but un- 
der the so-called strong-coupling conditions meaning in our notation y = 0.445 and x = 1.9 
{x/y = flji/uj2i = 4.27). Under these conditions our iterative spectrum generated from Eq. 
(|14p for = 30 cycle pulse is shown in Fig. (j2b)- An interesting feature of the iterative 
spectrum are (besides the familiar odd-order harmonics 3, 5 and 7) the doublets around the 
positions of even harmonics caused by the large Rabi splittings of the odd harmonics. This 
spectrum is a counter par t of the numerical spectrum of Zuo et al. (Fig. 7 in jl3|). 

Also Ivanov et al. |l6l [l7l | have calculated the emission spectra from molecular ions but 
using their analytical formula (Eq. (52) in jl^) derived in the limiting case of extremely 
strong coupling {y « 1 and x >> y in our language). We applied our Eq. (fT^ to this 
region and obtained with it the results shown in Fig. Q. This figure presents the heights 
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of the odd-harmonic peaks, H{n), normahzed to the third harmonic peak, for y = 0.1 and 
two values of x = 14.5 and 15, respectively. We have assumed a 30-cycle pulse of the form 
/(r)cosr with /(r) = 1. Our Fig. Q, obtained along the iterative procedure, coincides 
perfectly with the appropriate results of Ivanov and Corkum (Fig. 3 in [11]). 

A different place where two-level approximation has appeared to be reliable is a symmetric 
double quantum well E3, Hoi extensively studied in the context of laser control of 
tunneling and symmetry breaking with strong short pulses. The latter effect results in 
the appearance of spectral peaks at the positions of even harmonics from the systems with 
inversion symmetry. For example, Levinson et al.(Fig. 2 in their paper give the spectra 
from the double-well structure obtained by integrating numerically the set of three Bloch 
equations for the /(r) cosr pulse with /(r) = 1. The frequency-strength parameters in their 
numerical calculations fall into the strong-coupling region {y = 0.625, x = 1.25 in one case 
(their Fig. 2a) and y = 0.589, x = 1.178 in the other case (their Fig. 2b)). For the above 
two sets of frequency-strength parameters we show in Fig. (@)) our iterative spectra resulting 
from Eq. (jl4|) for the 30-cycle pulse assumed. The asymmetric doublets at the positions of 
the second and fourth harmonics, formed when using the first set of parameters (Fig. (^t-)); 
are seen to coalesce into single peaks when taking the other set (Fig. (jib)). Moreover, the 
second set of parameters results in shifting the low-frequency component of the spectrum 
towards zero. Both behaviours of our iteratively obtained spectra are the same as those in 
the numerical spectra of Levinson et al. (Fig. 2 in ^^) and are connected with approaching 
the so-called accidental degeneracy of two Floquet states of the system [l7, 3 2^ at some 
parameters. The parameter x from the second set does nearly satisfy the condition of the 
accidental degeneracy, i.e., it ensures for the Bessel function Jo(2x) to drop to zero [iTlli^. 

Miicke et al., have also used the two-level model to simulate numerically the spectra 
of light emitted around the third harmonic from GaAs semiconductor exposed to 5fs pulse 
of sech(r/ro) cosr form, where tq = tfwhm/I-^Q'^- The results of their simulations (Fig. 
2 in Q) reveal the evolution of the third-harmonic peak into a doublet structure when 
increasing the envelope pulse area A. For the sech(r/ro) envelope, the area A is related to 
our X parameter through A = titqx = {271^ / 1 .763) Npw h mx , where Nfwhm is the full width 
at half maximum (FWHM) measured in optical cycles {Nfwhm = 1-71 in this case). For 
the conditions close to those of Miicke et al., we present in Fig. (0) our iterative spectra, 
obtained from Eq. (fT^ . Our spectra are a qualitative reproduction of the numerical spectra 
of Miicke et al. (Fig. 2 in [^). A possible source of only qualitative agreement in this case 
is that our spectra are the pure response of the system, i.e., with no propagation effects 
included which were naturally taken into account in the simulations of Miicke et al. by 
coupling the Bloch equations to the Maxwell equations. 



IV. CARRIER ENVELOPE PHASE EFFECTS 

We now apply the iterative method to calculate, for a particular case, the dependence 
of the two-level-system response on the phase difference {(p ) between carrier wave and the 
maximum of the pulse envelope. To be specific, we make recurrence to the carrier-wave 
Rabi flopping originally studied numerically by Hughes 0] for a resonant {y = 1) pulse 
/i(r) = sech(r/ro) sin(r) of a fixed FWHM {Nfwhm = 1-72) but different pulse envelopes 
A = 19.24x. By coupling the optical Bloch equations to the Maxwell equations, Hughes 
considered propagation of the A = 2ln pulses through a two-level medium, where / was an 
integer. For the areas A = 6tt — 14tt of Hughes, the left-hand side column of Fig. [HI shows 
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population differences, — \b\l, versus time obtained by our iterative procedure witfi tfie 
use of Eq. fjll|) . Our results practically do not differ by nothing from the original numerical 
results obtained by Hughes just near the front-face of the two-level material, i.e^ where 
the propagation effects were not important yet (the left-side column of Fig. 3 in jj]). Our 
graphs confirm the original result of Hughes on incomplete Rabi flops at A > Svr. On the 
other hand, the right-side column of our Fig. IHl shows our iteratively obtained population 
differences but for the /i(r) = sech(r/ro) cos(r) pulse, i.e., the pulse with its carrier-envelope 
phase (p shifted by 7i/2 with respect to the pulse used by Hughes. Some differences caused by 
this shift are clearly seen in the middle parts of the population difference curves. These parts 
correspond to the times for which the pulse intensity has already evolved to its high values. 
The main differences introduced by changing the carrier-envelope phase (p consist in either 
converting the double peaks into single ones (and vice versa) or inverting the asymmetry in 
double peaks. 

The above 0-sensitivity of population inversion produces the dependence of the spectrum 
of scattered light on carrier-envelope phase. In Fig. 7a, we show the spectrum calculated 
iteratively with the use of Eq. (fT^ for the Hughes pulse /i(r) = sech(r/ro) sin(r + 0), i.e., 
with = 0, NpwHM = 1.72, y = I and x = 1.31 (this x corresponds to the envelope pulse 
area A = 8tt). Except the spectral peak at the fundamental frequency {z = 1), one sees 
a well pronounced peak at the position of second harmonic {z = 2) because the chosen 
X is in the vicinity of the value (1.178) ensuring the accidental degeneracy (Jo(2x) = 0) 
of the Floquet states of the system (compare Fig. HI and its discussion) . A similar peak 
around z = 2 was found by Tritschler et al. (Fig. la in jl^]) on the basis of their numerical 
solution of the two-level Bloch equations for a different pulse envelope (sinc(r/ro) instead 
of sech(r/ro)) and different light-matter parameters {NpwHM = 1-81, y = 2, x = 0.76). In 
addition to Fig. 7a, we show in Fig. 7b the dependence (calculated iteratively from Eq. 
(fT^) of the height of the spectral peak at 2; = 2 on the carrier-envelope phase < < 27r 
in the pulse h^r) = sech(r/ro) sin(r + 0). The 0-dependence is well seen and has a period 
of 71 in agreement with fully numerical calculations of Tritschler et al. 3| and Miicke et al. 



(Fig. lb in Igl) exploiting the optical Bloch equations. The same periodicity is seen in Fig. 
7c corresponding to the small peak in Fig. 7a around z = 1.5. 



V. SUMMARY 

On the basis of a nonlinear Riccati-type equation, analytically solved in very short time 
intervals (shorter than half the optical period), we have formulated an effective iterative pro- 
cedure for the problem of a two-level system exposed to a linearly polarized electromagnetic 
pulse. For different light-matter couplings (from weak through moderate to strong ones), 
we have proved very good quality of the procedure by recovering with it a number of subtle 
effects met in the previous numerically calculated photon-emission spectra and population 
inversion. We have applied the procedure developed to describe some carrier-envelope phase 
effects in extreme nonlinear optics, particularly in population inversion and spectrum of 
coherently scattered light. In the reg ime of a fe w-cy cle pulses, these carrier-envelope phase 
effects are of current interest H Q 111, H Q . 



If necessary, one could reinterpret the interaction Hamiltonian appropriately thus making 
the iterative procedure applicable to any form of time-dependent coupling within a two-level 
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system. 
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V SUMMARY 




Figure 1: The evolution of population inversion w versus K under the box-shaped (/(r) = 1) 
sine-like ((/> = 0) two-cycle {N = 2) pulse of x = y = 1. Solid lines - the results of the iterative 
procedure exploiting Eq. (Illj) . dotted lines - the results of direct numerical integration of the 
Riccati-type Eq. for r, dashed line - the results of square wave approximation. 
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Figure 2: Photon-emission spectra, from a model two-level molecular hydrogen ion, cal- 

culated iteratively by using Eq. dJ. (a): near-resonance case {y = 1.1 and x = 1.86) for the 
/(r) cos(t) pulse with /(r) gaussian increasing by 10 optical cycles and then keeping the value of 
1 up to 30 cycles; (b): strong-coupling case {y = 0.445 and x = 1.9) for the 30-cycle f (r) cosfr) 
pulse with / (t) = 1. The conditions are as those in the papers by Zuo et al. (Fig. 6b in [l5| and 
Fig. 7 in J^], respectively). For the high-resolution of the presented spectra we choose K = 100. 
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Figure 3: Normalized heights, H{n), of the odd-harmonic peaks calculated iteratively with the use 
of Eq. ((TH) . for the case of extremely strong coupling (y << 1, x >> y) of the two-level system 
to 30-cycle pulse of /(T)cos(r) form with /(t) = 1. The conditions are as those in the paper by 
Ivanov and Corkum (Fig. 3 in [l^l. 
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Figure 4: Spectra from a symmetric double-well structure obtained iteratively by applying Eq. 
((Til) , for 30-cycle pulse of /(r)cos(r) type with /(r) = 1 and frequency-strength parameters 
y = 0.625, X = 1.25 {(a)) and y = 0.589, x = 1.178 {(h)). The conditions are as those in the paper 
by Levinson et al. (Fig. 2 in 19]). 
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Figure 5: Evolution of the third-harmonic peak into a doublet with increasing pulse strength 
parameter x, calculated iteratively with the use of Eq. H15|). The conditions (5 fs pulse of 
sech (t/tq) cos (r) form) are close to those of numerical simulations by Miicke et al. (Fig. 2 in 
0). The positions of peaks in the doublet agree with those from numerical simulations but the 
heights are reversed. 
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Figure 6: Population difference, — |62p5 between the ground and excited states versus time, 
for /i(t) = sech(r/To) sin(r) pulse (left-side column) and /i(r) = secli(T/ro) cos(r) pulse (right-side 
column), i.e., the pulses with their carrier-envelope phases shifted to each other by 7r/2. The 
graphs were obtained iteratively by using Eq. (fTT|) . The chosen parameters: NpwHM = 1-72, 
y = 1, while the envelope pulse areas A and the corresponding strength parameters x are shown 
in the graphs. The left -column graphs reproduce the fully numerical results of Hughes (Fig. 3 in 
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Figure 7: (a) The iteratively calculated (from Eq. (jlSI) ) photon emission spectrum generated by 
/i(r) = sech(r/ro) sin(T + cp) pulse of the parameters </> = 0, NpwHM = 1-72, y = 1 and x = 1.31 
(Stt pulse). (6) The calculated (from Eq. ^(T^ ) height of the spectral peak at the position of second 
harmonic (z = 2) versus carrier-envelope offset phase 4>. The (^-dependence with a period of vr 
agrees with the results of fully numerical calculations by other Authors 0, for different pulse 
shapes and light-matter parameters, (c) The same as (b) but for the small peak at z = 1.5 in (a). 



